Citation: CPT Pharmacometrics Syst. Pharmacol. (2014) 3, el03; doi:10.1038/psp.2013.80 
© 2014 ASCPT All rights reserved 2163-8306/14 



www.nature.com/psp 

ORIGINAL ARTICLE 

Population Pharmacokinetics of Azithromycin in 
Whole Blood, Peripheral Blood Mononuclear Cells, 
and Polymorphonuclear Cells in Healthy Adults 

MR Sampson 1 , TP Dumitrescu 12 , KLR Brouwer 1 and VD Schmith 2 

Azithromycin's extensive distribution to proinflammatory cells, including peripheral blood mononuclear cells (PBMCs) 
and polymorphonuclear cells (PMNs), may be important to its antimicrobial and anti-inflammatory properties. The need to 
simultaneously predict azithromycin concentrations in whole blood ("blood"), PBMCs, and PMNs motivated this investigation. 
A single-dose study in 20 healthy adults was conducted, and nonlinear mixed effects modeling was used to simultaneously 
describe azithromycin concentrations in blood, PBMCs, and PMNs (simultaneous PK model). Data were well described by a 
four-compartment mamillary model. Apparent central clearance and volume of distribution estimates were 67.3 l/hour and 336 1 
(interindividual variability of 114 and 122%, respectively). Bootstrapping and visual predictive checks showed adequate model 
performance. Azithromycin concentrations in blood, PBMCs, and PMNs from external studies of healthy adults and cystic 
fibrosis patients were within the 5th and 95th percentiles of model simulations. This novel empirical model can be used to 
predict azithromycin concentrations in blood, PBMCs, and PMNs with different dosing regimens. 
CPT Pharmacometrics Syst. Pharmacol. (2014) 3, e103; doi:1 0.1 038/psp.201 3.80; published online 5 March 2014 



Azithromycin is approved worldwide as a broad-spectrum 
antibiotic to treat a variety of community-acquired infec- 
tions. The recommended dose of azithromycin for the treat- 
ment of acute respiratory tract infections is 500 mg (10mg/ 
kg in children) once daily for 3 days, or 500 mg (10mg/ 
kg) on day 1, followed by 250mg (5mg/kg) once daily for 
4 days. 1 The major route of elimination is biliary excretion. 
Azithromycin and other macrolides are under investigation 
for anti-inflammatory/immunomodulatory actions that may be 
beneficial in chronic neutrophil-driven pulmonary inflamma- 
tory syndromes. 2-4 

The mechanism of anti-inflammatory action of macrolides 
is unknown. Macrolides affect a variety of PMN (or neu- 
trophil) functions in vivo, including chemotaxis, adhesion, 
degranulation, cytokine release, and apoptosis. 56 In addi- 
tion, activated mononuclear phagocytes (a subset of periph- 
eral blood mononuclear cells (PBMCs), which also include 
lymphocytes and macrophages) may be cellular targets of 
their anti-inflammatory action. 7 Interestingly, there is no clear 
in vitro/in vivo correlation, and short-term direct effects on 
PMNs cannot explain the fact that in chronic inflammatory 
conditions, weeks to months of macrolide therapy appears 
to be required for symptomatic improvement. 8 For example, 
a pharmacokinetic (PK)/pharmacodynamic (PD) model pre- 
dicting inhibition of interleukin-6 and tumor necrosis factor-a 
after high intravenous doses of erythromycin 9 predicts little or 
no inhibition of interleukin-6 and tumor necrosis factor-a at 
250 mg once daily, a dose that produces reductions in exac- 
erbations of bronchiectasis. 10 

The optimal dosing regimen for chronic treatment with 
azithromycin is unknown. Numerous empiric dosing regi- 
mens (once daily, every other day, once weekly, twice weekly, 



thrice weekly) have been studied to evaluate the anti-inflam- 
matory effects of azithromycin in asthma, chronic obstructive 
pulmonary disease, cystic fibrosis (CF), and bronchiecta- 
sis. 11-15 Differences in dosing regimens are unlikely to be 
explained by PK differences between populations given that 
azithromycin PK between healthy volunteers and CF patients 
were similar. 16 

Azithromycin is taken up quickly into human fibroblasts, 
PBMCs, and PMNs, stored in lysosomes, and slowly released, 
most likely by a pH-dependent process. 1718 Thus, intracellu- 
lar retention of azithromycin is expected, with azithromycin 
concentrations in monocytes (a subtype of PBMCs) still pres- 
ent 14 days after 5 days of dosing in healthy volunteers. 17 
Higher concentrations of azithromycin have been observed 
in infected vs. uninfected tissue in a mouse thigh infection 
model and in inflamed vs. noninflamed blisters in humans. 1920 
The contribution of intracellular vs. extracellular macrolide 
concentrations to efficacy is also unknown. At a minimum, 
the intracellular uptake, retention, and release of azithro- 
mycin may improve drug delivery to the site of infection or 
inflammation. 6 ' 1819 ' 21 Therefore, an understanding of the rela- 
tionship between azithromycin concentrations in whole blood 
(hereafter referred to as "blood"), PBMCs, and PMNs may be 
more relevant to efficacy than plasma concentrations. Future 
studies aimed at understanding the relationship between 
intracellular and blood azithromycin concentrations and effi- 
cacy would be ideal. However, considering the poorly under- 
stood mechanism of action for azithromycin and the lack of 
PK/PD relationships, the most practical route to optimizing 
azithromycin intermittent dosing will be based on a well- 
evaluated PK model of azithromycin concentrations in blood, 
PBMCs, and PMNs over time. 
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Several PK studies have reported azithromycin concentra- 
tions in plasma, serum, or blood in addition to PBMCs and/ 
or PMNs in volunteers, 202223 CF patients, 162425 and human 
immunodeficiency virus-infected patients. 26 Limitations of 
these studies included a short sampling collection inter- 
val (e.g., 3-7 days) relative to the reported half-life (up to 
200 hours in blood) 25 and the fact that these studies did not 
include models characterizing simultaneous azithromycin 
concentrations in blood, PBMCs, and PMNs over time. 

Blood is a mixture of red blood cells, which account for 
40-50% of the total volume, white blood cells (including 
PBMCs and PMNs), which account for -1% of the total vol- 
ume, and plasma, which accounts for the remaining 55% of 
the blood volume. Of the white blood cells, PBMCs account 
for 28-40%, and PMNs account for 54-62%. 27 Azithromycin 
concentrations are relatively low in plasma and blood, mini- 
mal in red blood cells, and >1 00-fold higher than blood in 
PBMCs and PMNs. Because the mechanism of anti-inflam- 
matory action of macrolides is unknown, there are different 
opinions as to whether PBMCs or PMNs are more important 
to macrolide efficacy. Thus, we evaluated both PBMCs and 
PMNs in this study. The present study was designed to: (i) 
characterize the PK of azithromycin in blood, PBMCs, and 
PMNs for 3 weeks after a single dose of azithromycin and (ii) 
develop a simultaneous PK model that could be used to pre- 
dict long-term human exposure to azithromycin and/or novel 
compounds with similar PK properties. 

RESULTS 

Study population and samples 

Twenty participants with a median age of 48.5 years (range: 
21-63) and body mass index of 24.9 kg/m 2 (21.4-28.2) 
were enrolled. Twelve (60%) participants were male and 16 
(80%) were Caucasian. A total of 269 blood, 227 PBMC, 
and 239 PMN concentrations were included in the analysis. 
One PBMC sample was excluded due to issues with sample 
processing. 
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Figure 1 Azithroycin mean (±SE) (a,b) and individual (c,d) 
concentration-time profiles for all timepoints (a,c) and for the first 
24-h period (b,d). Dashed line, 250-mg dose group; solid line, 
1 ,000-mg dose group; green, whole blood; red, peripheral blood 
mononuclear cells; blue, polymorphonuclear cells. 



Concentration-time profiles for each participant are shown 
in Figure 1 , stratified by dose. Azithromycin blood concentra- 
tions peaked early (median T max of 3.5 hours) and declined 
biexponentially with time, whereas azithromycin concentra- 
tions in PBMCs and PMNs were much more variable, with a 
later peak (median T max of 9.0 hours in PBMCs and PMNs), 
as well as multiple, inconsistent, unsystematic fluctuations in 
individuals, particularly after 48 hours (Figure 1). For the 10 
participants receiving the 1,000-mg dose, azithromycin con- 
centrations in blood, PBMCs, and PMNs were measurable in 
6, 9, and 7 participants (median concentrations of 0.01, 7.14, 
and 2.18mg/l) at the 3-week timepoint, respectively. 

Population PK model development 

Numerous models were explored to describe the fluctuations 
in azithromycin concentrations in PBMCs and PMNs. An 
oscillating sine function (see Supplementary Equation S1) 
and modified enterohepatic recirculation model was applied 28 
to attempt to characterize fluctuations due to the rapid uptake 
of azithromycin into, and very slow release from PBMCs and 
PMNs, with further reuptake. As an attempt to account for 
uncertainty in cell counts or cell volumes, cell count and 
cell volume were included in the residual error model. None 
of these approaches resulted in adequate goodness of fit. 
Therefore, empirical strategies were evaluated. 

The final model that best described all of the data was a 
four-compartment, mamillary, empirical model with first-order 
absorption, first-order elimination from the central compart- 
ment (Compl), unidirectional flow between the central com- 
partment and cellular compartments (Comp2 for PBMCs 
and Comp3 for PMNs), bidirectional distribution between the 
central compartment and large peripheral compartment ("tis- 
sue/rest of body", Comp4), and elimination from Comp1-3 
(Figure 2; see Supplementary Materials for model code). 
Unidirectional flow between the central and cellular compart- 
ments was used because a fixed azithromycin uptake-to- 
efflux ratio (based on in w'frodata in PMNs) 17 with bidirectional 
flow (where cellular uptake rate » cellular release rate) did 
not provide adequate goodness of fit. 
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Figure 2 Four-compartment mamillary structural model. Compl , 
central compartment; Comp2, peripheral blood mononuclear cell 
compartment; Comp3, polymorphonuclear cell compartment; and 
Comp4, tissue compartment. Azithromycin concentrations in whole 
blood are equal to A*[Comp1] + S*[Comp2] + C[Comp3], where the 
coefficients A and 6 were estimated empirically, and Cwas fixed. 
CL12, 13, 14, and 41 are intercompartmental clearances; CL1-3 
are elimination clearances. 
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As blood contains plasma, PBMCs, and PMNs, azithro- 
mycin concentrations in blood were modeled as the sum of 
concentrations in each compartment multiplied by a coeffi- 
cient representing the fraction of blood occupied by plasma, 
PBMCs, and PMNs, respectively: [Blood] = /4[Comp1] + 
6[Comp2] + C[Comp3]. After unsuccessfully attempting to 
include physiologically relevant parameter estimates for the 
coefficients, A and B were estimated empirically, and C was 
fixed for the final model. Proportional residual variability and 
exponential interindividual variability models were used. The 
data did not support the addition of any covariates. 

Model evaluation 

In each matrix, approximately symmetrical distributions of 
observed vs. population (PRED) and individual (IPRED) pre- 
dicted concentrations about the line of unity, and lack of bias 
in conditional weighted residuals (CWRES) vs. PRED and 
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Figure 3 (a) Population predicted and (b) individual predicted vs. 
observed azithromycin concentrations. Green circles, whole blood; 
red diamonds, peripheral blood mononuclear cells; blue triangles, 
polymorphonuclear cells. 
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Figure 4 Conditional weighted residuals vs. (a) population 
predictions (PRED) and (b-d) time. Green circles, whole blood; 
red diamonds, peripheral blood mononuclear cells; blue triangles, 
polymorphonuclear cells; dashed line, lowess trend line. 



time, indicated adequate goodness of fit (Figures 3 and 4). 
The median (range) percent difference between the model 
and bootstrap median estimates was 5% (0-18), and rela- 
tive standard errors of the fixed effect parameters were all 
<20%, indicating good precision of parameter estimates 
(Table 1). With the exception of the interindividual variability 
on the absorption rate constant (shrinkage = 52%), shrink- 
ages were <30%. The random effect on the absorption rate 
constant was retained in the model to better characterize the 
azithromycin blood concentrations at T max . In the visual pre- 
dictive checks (VPCs) for blood, PBMCs", and PMNs, 12.3, 
9.7, and 6.7% of observed concentrations were outside the 
90% prediction interval, respectively, indicating a good over- 
lap between the observed and simulated concentrations (see 
Figure 5 for data up to 3 weeks and Supplementary Figure 
S1 for first 24 hours). External data from healthy volunteers 
and CF patients were within the 5th and 95th percentiles of 
model simulations (Figure 6). 

Normalized prediction distribution errors (NPDEs) vs. pre- 
dicted concentrations (YPRED) and vs. time plots were eval- 
uated post hoc. These plots showed bias early in the NPDE 
vs. YPRED plots and late in the NPDE vs. time plots for all 
matrices. As NPDEs are a relatively new method that takes 
into account the correlation in residuals within individuals, it 
is unknown how this method performs when applied to data 
sets with high intrinsic measurement variability in small num- 
bers of participants, such as the data set in the present study. 
The bias observed in the NPDE plots are not supported by 
the CWRES, VPC, and bootstrapping results, and therefore 
are unlikely to significantly impact model-based decision 
making. 



DISCUSSION 

Azithromycin distribution between blood, PBMCs, and PMNs 
was well described by a four-compartment model. Internal 
and external evaluation indicated good model performance 
for single- and multiple-dose regimens and good model pre- 
dictions compared with reported external data in healthy 
adults and CF patients. Thus, the present simultaneous PK 
model can be used to predict the long-term exposure in 
humans receiving various azithromycin intermittent dosing 
regimens. In addition, this model can be used for animal-to- 
human translation of relevant dosing regimens (e.g., once 
daily, once weekly) for novel compounds exhibiting azithro- 
mycin-like PK. 

Multiple, inconsistent, unsystematic fluctuations in azithro- 
mycin concentrations in PBMCs and PMNs over time were 
observed in individual participants, particularly after 48 hours. 
Azithromycin is known to be sequestered in PBMCs, PMNs, 
and fibroblasts, 1829 and the average life span of PBMCs 
and PMNs are a few weeks and a few days, respectively. 30 
Thus, we initially hypothesized that the observed fluctuations 
resulted from azithromycin release due to cell turnover, fol- 
lowed by reuptake by other cells. However, it is unlikely that 
all PBMCs and PMNs would release and reuptake azithro- 
mycin in a coordinated manner and result in the fluctuations 
observed. Attempts to account for these fluctuations with an 
oscillating function or as an artifact of sample processing 
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Table 1 Model and bootstrap parameter estimates for the final model 

Bootstrap 
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CL13/F(l/h) 
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31.4 


75.4 


CL3/F(l/h) 


0.026 
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Proportional residual variability and exponential random effect models 
were used. rjCL1 2-CL1 3/F is the shared ti estimate for CL12 and CL13; 
C=B/1,000; CL41 =CL14/2. 

PBMC, peripheral blood mononuclear cell; PMN, polymorphonuclear cell 
RSE, relative standard error. 
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Figure 5 Visual predictive checks. Markers, observed data; shaded 
area, 90% prediction interval; solid lines, median of observed data; 
dashed lines, median of simulated data; green circles, whole blood; 
red diamonds, peripheral blood mononuclear cells; blue triangles, 
polymorphonuclear cells. 



in the residual variability model did not improve goodness 
of fit. Therefore, the final model presented herein assumed 
that these fluctuations were due to random variability. Future 
studies are needed to rule out whether fluctuating azithromy- 
cin concentrations in PBMCs and PMNs result from a biologi- 
cal mechanism. 

The magnitudes of interindividual variability on V1/F, V2/F, 
1/3/F, CL1/F, and CL12-CL1 3/F were relatively large. The 
small and relatively homogenous patient population did not 
support inclusion of covariates. Nevertheless, our model pre- 
dictions for azithromycin concentrations in PBMCs and PMNs 
were in agreement with concentration-time data reported in 
prior studies of healthy volunteers and CF patients. The com- 
parable concentrations observed between CF patients and 
predictions from our healthy volunteer model are consistent 
with a study comparing azithromycin PK in 12 healthy volun- 
teers and 12 patients with CF, which reported no PK differ- 
ences between groups. 16 

One previous oral azithromycin PK study in healthy adults 
measured azithromycin blood concentrations and used a popu- 
lation modeling approach. 31 Using a three-compartment model, 
central CL and total V(sum of all V terms) were 115 l/hour and 
6,238 I, respectively, compared with the estimates in the pres- 
ent study (four-compartment model) of 67.3 l/hour and 4,937 
I, respectively. The differences in CL and V estimates between 
studies may be related to differences in sampling times. The 
present study had more frequent sampling times over a 3-week 
period, whereas the previous study had multiple samples over 
the first 48 hours and then no samples until 3 weeks. 

For over two decades, azithromycin's extensive intracellu- 
lar and tissue distribution have been known, and proposed to 
be related to efficacy. 1832 However, target azithromycin con- 
centrations in inflammatory cells (i.e., PBMCs and PMNs) for 
anti-infective or anti-inflammatory uses have not been estab- 
lished. For example, systemic and intracellular azithromycin 
PK has been evaluated in small studies in CF, 1624 25 but larger 
efficacy studies in CF have not collected PK information that 
would enable evaluation of target concentrations. 33-35 

Blood is a mixture of red blood cells, plasma, and white 
blood cells. As there is minimal uptake of azithromycin into red 
blood cells, and extensive uptake and retention of azithromy- 
cin in PBMCs and PMNs, 36 the interpretation of blood azithro- 
mycin concentrations over time is complex. It is important to 
note that at early timepoints, blood concentrations parallel 
plasma concentrations 31 ; however, during the terminal phase, 
concentrations of azithromycin in blood parallel concentra- 
tions of azithromycin in PBMCs and PMNs (which comprise a 
small fraction of total blood volume, but contain much higher 
azithromycin concentrations). Since isolation of PBMCs and 
PMNs are labor intensive, simulations from this model could 
provide information on the best timing for limited sampling of 
PBMCs and PMNs in future studies. 

In conclusion, a model that simultaneously described 
azithromycin concentrations in blood, PBMCs, and PMNs 
over time was developed with good model performance. This 
model can be used to predict azithromycin concentrations in 
blood, PBMCs, and PMNs following various oral dosing regi- 
mens and to aid in selection of the most appropriate dosing 
regimens for chronic administration of azithromycin or novel 
compounds with azithromycin-like PK. 
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Figure 6 Comparison of model predictions to external data. 
Azithromycin concentrations in (a) peripheral blood mononuclear 
cells (PBMCs) and (b) polymorphonuclear cells (PMNs) of healthy 
individuals after receiving 500 mg of azithromycin daily for 3 days. 
Small solid markers with error bars, reported azithromycin PBMC 
concentrations in 24 healthy volunteers (mean ± SD); large 
open markers without error bars, reported mean azithromycin 
concentrations in 12 healthy volunteers; red diamonds, PBMCs; 
blue triangles, PMNs; solid lines represent the 5th, 50th, and 95th 
percentiles of model-simulated (n = 5,000) concentrations, (c) C max 
and (d) C min of azithromycin concentrations in PMNs in cystic fibrosis 
(CF) patients after receiving 500 mg of azithromycin daily for 35 days. 
Open blue triangles, reported azithromycin PMN concentrations in 
eight CF patients (mean ± SD); solid lines represent the 5th, 50th, 
and 95th percentiles of model-simulated (n = 5,000) concentrations. 

METHODS 
Study design 

A prospective, single-dose, PK trial in healthy male and 
female volunteers was conducted. Inclusion criteria included 
the following: healthy based on medical history, physical 
examination, laboratory tests, and cardiac monitoring; age 
of 18-65 years; willing to use contraception throughout the 
study (women of child-bearing potential); body weight >50kg 
and body mass index within the range of 1 8.5-30 kg/m 2 ; 
provided informed consent; and corrected QT interval 
(electrocardiogram time interval from the start of the Q wave 
to the end of the T wave) <450 ms. Exclusion criteria included 
the following: current or history of positive prestudy hepati- 
tis B surface antigen or positive hepatitis C antibody result; 
current or chronic history of liver disease or known hepatic 
or biliary abnormalities (with the exception of Gilbert's syn- 
drome or asymptomatic gallstones); a positive prestudy 
drug/alcohol screen; a positive test for human immunodefi- 
ciency virus antibody; history of regular alcohol consumption 
within 6 months of the study; participation in an investiga- 
tional product trial within the last 3 months; exposure to more 
than four new chemical entities within the last 1 2 months; use 
of prescription or non-prescription drugs, including vitamins, 
herbal, and dietary supplements within 7 days (or 14 days 
if the drug is a potential enzyme inducer) or five half-lives 
(whichever is longer); history of sensitivity to azithromycin or 
other macrolides; donation of blood in the past 3 months; 
pregnant or lactating females; unwillingness or inability to fol- 
low the procedures outlined in the protocol; and mentally or 
legally incapacitated. 



Participants were randomized to receive a single oral 
azithromycin dose of 250 or 1,000mg. Serial blood samples 
were collected predose, and at 1, 2, 3, 4, 6, 9, 12, 16, 24, 48, 
96, 144, 240, 336, and 504 hours postdose. Azithromycin con- 
centrations in blood were measured at all timepoints, whereas 
concentrations in PBMCs and PMNs were measured at all 
timepoints except the 3- and 240-h timepoints. Blood volumes 
of 0.5 and 10 ml per timepoint were collected for azithromy- 
cin concentration determination in blood and PBMCs/PMNs, 
respectively. The study was approved by the Surrey Research 
Ethics Committee, conducted in accordance with the ethical 
principles set forth in the Declaration of Helsinki, and posted to 
the clinicaltrials.gov (NCT01416350) website. 

Sample processing 

PBMCs and PMNs were isolated from hemolyzed whole 
blood using density gradient centrifugation. Five milliliters 
of Histopaque-1 1 19 (Sigma Aldrich, St Louis, MO) at room 
temperature was layered slowly into a tube containing 5 ml 
room temperature Histopaque-1 077 (Sigma Aldrich). Ten mil- 
liliters of anticoagulated blood was slowly layered onto the 
Histopaque tube and the tube was centrifuged at 700gfor 30 
minutes at 20 °C with the brake off. The PBMCs in the opaque 
layer just above the Histopaque-1 077 were transferred to a 
fresh tube. The PMNs in the opaque layer just above the His- 
topaque-1 119 were transferred to an additional fresh tube. 
At this point, each PBMC and PMN sample was processed 
identically. Each tube was filled to 12 ml with cold phosphate- 
buffered saline and centrifuged at 200g> for 10 minutes at 4 °C 
with the brake on. The supernatant was discarded, and the 
cells were promptly resuspended. Tubes were filled to 12 ml 
with phosphate-buffered saline and centrifuged at 200g for 
10 minutes at 4 °C with the brake on. The supernatant was 
discarded and the pellet was resuspended in exactly 10 ml of 
phosphate-buffered saline and thoroughly mixed. A 200 pi of 
aliquot was taken to measure the cell count using a Horiba 
ABX P60 (Horiba, Irvine, CA) hematology analyzer. The sam- 
ple was aliquoted into fresh tubes containing <1 x10 7 cells 
per tube. These tubes were filled to 12 ml with phosphate- 
buffered saline and centrifuged at 200g for 10 minutes at 4 °C 
with the brake on. The supernatant was decanted and 150 pi 
of Triton lysis buffer was added. Samples were mixed by vor- 
tex, weighed to determine the final volume of the cell lysate, 
and then frozen at -80 °C. 

Bioanalytical assay 

Azithromycin concentrations in hemolyzed whole blood, 
and PBMC and PMN lysate samples were measured using 
a validated tandem mass spectrometry method for concen- 
trations between 0.002 and 2 mg/l. A Waters Acquity UPLC 
(Waters, Milford, MA) with a mobile phase flow rate of 
0.8 ml/minute and an Applied Biosystems/MDS Sciex API- 
4000 mass spectrometer (Applied Biosystems, Foster City, 
CA) with a TurbolonSpray (Applied Biosystems) interface 
and multiple reaction monitoring were used. Azithromycin 
was extracted from 25 pi of human blood by protein precipi- 
tation using acetonitrile containing an isotopically labeled 
internal standard ([ 2 H 5 ]-azithromycin). The same method 
also was validated for the determination of azithromycin 
in PBMC and PMN lysate samples, using human blood 
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calibration standards. In validation runs for blood, the maxi- 
mum bias was 9.7% and the maximum within-run preci- 
sion value was 14.4% at concentrations of 0.002, 0.006, 
0.1, 1.6, and 2mg/l. In PBMC and PMN validation runs, 
the maximum bias was 14.9% and the maximum within-run 
precision value was 8.1% at concentrations of 0.006 and 
1.6mg/l. As the bias and precision values were all <15% 
(<20% at lower limit of quantification), the method was in 
accordance with bioanalytical guidelines. 

Calculation of intracellular azithromycin concentrations 

Azithromycin concentrations in PBMCs and PMNs were cal- 
culated using the following equations, assuming an average 
PBMC cell volume (0.204 pi) and PMN cell volume (0.334 pi), 
as reported previously 37 : 

Total volume of cells (ul) = number of cells in pellet 
x10 6 x average cell volume (pi) x1 ul/10 6 pl 

Azithromycin intracellular concentration (mg/l) 

= azithromycin concentration in lysate(mg/l) (2) 

x lysate volume (ul) / total cell volume (ul) 

Population PK model development 

Azithromycin concentration-time data in blood, PBMCs, and 
PMNs were analyzed simultaneously using nonlinear mixed 
effects modeling with Phoenix NLME 2.1 (Certara, St Louis, 
MO) software, utilizing the first-order conditional estimation 
method with interaction. Various structural compartmental 
mamillary and catenary models with uni- or bidirectional 
intercompartmental transfer, coupled with various random 
effects models, were explored. 

Within-participant fluctuations in azithromycin concentra- 
tions in PBMCs and PMNs were first modeled with a sine 
function, which has been used to describe periodic increases 
in concentration due to reabsorption (e.g., for a drug that 
undergoes enterohepatic recirculation). 28 This equation (see 
Supplementary Equation S1), intended to describe fluctu- 
ating, or periodic change in the drug transfer rate from one 
compartment to another (e.g., drug leaving PMNs, returning 
to plasma or some tissue compartment, and then being taken 
up again into PMNs), was then incorporated into one or more 
structural model differential equations. 

Another approach used to account for fluctuations in 
azithromycin concentrations in PBMCs and PMNs was to 
incorporate PBMC and PMN cell volume or cell count (mea- 
sured during sample processing) into the residual error 
model. This residual error model was chosen to account 
for the complex nature of processing these cells (separa- 
tion of cells, counting of cells, and determination of lysate 
concentration) and the potential for inherent variability in 
the process. 

The appropriate model was chosen based on good- 
ness of fit plots, plausibility of parameter estimates, Akaike 
information criterion, a reduction in residual variability, and 



precision of parameters. Diagnostic plots included the follow- 
ing: observed concentrations vs. PRED and IPRED; CWRES 
vs. PRED and vs. time after dose; CWRES histogram; histo- 
grams of each random effect; and observed vs. PRED over 
time by individual. Random effects were considered to be 
supported by the data if shrinkage was <30%. 

Model evaluation 

Candidate final models were bootstrapped (n = 1,000) to 
assess the precision of parameter estimates. The relative dif- 
ferences between final model and bootstrap parameter esti- 
mates were calculated as follows: 

Percentage difference = (F-S)/Fx100% (3) 

where Fis the final model parameter estimate and B is the 
bootstrap parameter estimate. VPCs (n = 1 ,000 simulations) 
were conducted for candidate final models to compare the 
degree of overlap between observed and model-simulated 
concentrations. External published azithromycin concen- 
trations in PBMCs or PMNs obtained from studies of 24 
healthy volunteers, 20 12 healthy volunteers, 25 and 8 CF 
patients 24 were digitized and used for external evaluation of 
the final model. These data were overlaid with the 5th, 50th, 
and 95th percentiles of the model-simulated concentra- 
tion-time curves for the dosing regimen used in each study 
(n = 1,000-5,000 simulations). Bootstrapping and simula- 
tions were conducted using Phoenix NLME 2.1. VPCs and 
simulations were plotted using Stata 12 (StataCorp, Col- 
lege Station, TX). Data were digitized using Data Digitizer 
(http://plotdigitizer.sourceforge.net/). NPDE vs. YPRED 
and NPDE vs. time plots were evaluated post hoc (follow- 
ing selection of the final model). The simulations used for 
VPCs also were used for computation of NPDEs. NPDEs 
and YPRED were calculated using an R software add-on 
package. 38 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y Intracellular concentrations of azithromycin in 
PBMCs and PMNs may be relevant to the an- 
timicrobial and anti-inflammatory properties of 
this drug. Most studies evaluating azithromycin 
intracellular concentrations have focused on 
only one cell type, collected sparse samples, 
or measured concentrations for only 3-7 days. 
Azithromycin has shown promising efficacy in 
chronic neutrophil-driven inflammatory diseas- 
es, although optimal dosing is unknown. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

y This study characterized the relationship be- 
tween azithromycin concentrations in whole 
blood, PBMCs, and PMNs for 3 weeks after a 
single dose. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

S This is the first PK study to intensively sample 
intracellular azithromycin concentrations and si- 
multaneously model the time course of azithro- 
mycin in systemic and intracellular matrices 
over a long (3-week) time course. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y The model enables assessment of azithromycin 
systemic and intracellular exposure under vari- 
ous dosing regimens for chronic investigational 
anti-inflammatory conditions and enables as- 
sessment of an optimal sparse sampling strat- 
egy for PBMCs and PMNs. 
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